Quantum many-body theory of qubit decoherence in a finite-size spin bath. II. Ensemble dynamics 
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Decoherence of a center spin or qubit in a spin batii is essentially determined by the many-body bath evolution. 
The bath dynamics can start either from a pure state or, more generally, from a statistical ensemble. In the 
preceding article [W. Yang and R. B. Liu, Phys. Rev. B 78, 085315 (2008)], we have developed the cluster- 
correlation expansion (CCE) theory for the so-called single-sample bath dynamics initiated from a factorizable 
pure state. Here we present the ensemble CCE theory, which is based on similar ideas of the single-sample CCE: 
The bath evolution is factorized into the product of all possible cluster correlations, each of which accounts for 
the authentic (non-factorizable) collective excitation of a group of bath spins, and for the finite-time evolution in 
the qubit decoherence problem, convergent results can be obtained by truncating the ensemble CCE by keeping 
cluster correlations up to a certain size. A difference between the ensemble CCE and single-sample CCE is that 
the mean-field treatment in the latter formalism of the diagonal part of the spin-spin interaction in the bath is not 
possible in the former case. The ensemble CCE can be applied to non-factorizable initial states. The ensemble 
CCE is checked against the exact solution of an XY spin bath model. For small spin baths, it is shown that 
single-sample dynamics is sensitive to the sampling of the initial state from a thermal ensemble and hence very 
different from the ensemble average. 

PACS numbers: 76.20.+q, 03.65.Yz, 76.60.Lz, 76.30.-v 



I. INTRODUCTION 

The dissipative dynamics of a center spin in a spin bath' 
is an old topic in spin resonance spectroscopy-^'^ '*'^ Recently, 

this subject is revisited^-^i^i^i'°-"i'^i'3.'4.'5,16i'7,lSi'9,20,21,22.23,24 

mostly due to the decoherence issue in quantum informa- 
tion processing. ^^-^^'^^ Being a most promising candidate for 
solid state qubits, electron spins in quantum dots or impu- 
rity centers experience decoherence by coupling to complex 
solid-state environments. A series of theoretica l^-"*'"''^ and 
experimental^'^ *^ works have identified that the dominating 
decoherence mechanism for electron spin qubits at low tem- 
peratures (such as below a few Kelvins) is the entanglement 
with nuclear spins of the host lattice i'^-''*''^-"''''^-'^-'V°'^' 

When the qubit flip is suppressed (usually by the large Zee- 
man energy mismatch between qubit and bath spins in a mod- 
erate magnetic field), the Hamiltonian for a qubit-bath system 
has the general form 

i/ = |-f)i/W(+| + |_>//(-)(-|. (1) 

The bath dynamics is driven by different Hamiltonians W'^^ 
depending on the qubit states |+). For a given initial bath state 
IJT) (which could be a random sampling from a thermal en- 
semble), the qubit coherence at time T is characterized by the 
"single-sample" propagator (^le''^' ''e^'^^*''\J'). For a thermal 
ensemble of bath states characterized by a density matrix p, 
a further ensemble average should be processed and the qubit 
coherence is given by the ensemble average Tr ype'^^ j 
In general, the key is to evaluate the ensemble-averaged prop- 
agator 

£ = Tr(p.'■«'^'■«'^..) (2) 

for a general density matrix p and arbitrary bath interaction 
operators {(9^-''). To address this problem, a variety of quantum 
many-body theories have been developed, including the den- 
sity matrix cluster expansion j'-^''"*''^-'^''^ the pair-correlation 



approximation ^ 1^' " and the linked-cluster expansionist The 
pair-correlation approximation provides a clear physical pic- 
ture for the bath dynamics by keeping only spin-pair corre- 
lations. The linked-cluster expansion accurately takes into 
account higher-order correlations with a Feynman diagram 
method, which, however, becomes dramatically tedious with 
increasing the order of diagrams. The density matrix cluster 
expansion simplifies the evaluation of higher-order correla- 
tions, but it may not converge to the exact results for relatively 
small bathsi^ 

Very recently, we have developed a cluster-correlation ex- 
pansion (CCE) theory^"* for the evaluation of the single- 
sample propagator Lj = (JTIe'''' e'''' ' ■ ■ ■ \ J'), which is a spe- 
cial case of the ensemble-averaged propagator in Eq. (|2]i. The 
CCE method provides a simple and accurate method to sys- 
tematically take into account the high-order correlations. For 
a temperature T much higher than the bath interaction strength 
(~10 K for nuclear spins in GaAs), the initial thermal en- 
semble can be well approximated as p a; expi-HQ/ksT), 
where Ho is the non-interacting Hamiltonian containing only 
the Zeeman energy. Such an initial ensemble is factorizable 
and a sampling \ from the ensemble can be taken as a prod- 
uct state \ J'} - Ijn) of all constituent bath spins, where 
Ijn) denotes the Zeeman energy eigen state of the nth bath 
spin. For a large spin bath, previous study^° has shown that 
the qubit decoherence is insensitive to the random sampling 
of the initial bath state from a thermal ensemble since the sta- 
tistical fluctuation scales with the number of bath spins as 
1 / Va. Thus the ensemble dynamics can be just identified 
with the single-sample dynamics with a random choice of the 
initial state. For a relatively small bath, however, the single- 
sample dynamics could be sensitive to the sampling of the ini- 
tial state and the ensemble average can be very different from 
any single sample. More importantly, if the initial state of the 
bath is entangled, i.e., 

P*T,P»®Ptv (3) 

a i 
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for any choice of probability distribution {Pa} and single spin 
density matrices {p|)^^), the single-sample CCE is not appli- 
cable. To extend to general ensemble bath dynamics, one 
could simply use the Monte Carlo simulation with a suffi- 
ciently large random sampling of the initial states from the 
ensemble. The Monte Carlo simulation is practically cumber- 
some due to the large number of initial states required for a 
faithful reproduction of the ensemble dynamics, and more im- 
portantly, it cannot be applied to non-factorizable initial states. 
In this paper we will develop a CCE formalism suitable for di- 
rect evaluation of ensemble-averaged bath evolution. 

In Sec. II, we will present the ensemble CCE and compare 
it with the single-sample CCE. In Sec. Ill, we check the en- 
semble CCE against the exact solution of a one-dimensional 
XY mode and compare the single-sample CCE and ensemble 
CCE. Sec. IV gives the conclusions. 



II. ENSEMBLE CLUSTER-CORRELATION EXPANSION 

A. An example 

Let us consider a bath consisting of spins and evaluate 
the bath evolution 

£sTr(pe'0), (4) 
averaged over a noninteracting (factorizable) ensemble 

p = P|l) (8)p|2| ® ■ ■ ■ ®P|Af), (5) 

where p(i) = 2; Pj\j)ij\ is the non-interacting density matrix 
for the rth spin J, and 

n m<n 

is the dimensionless bath interaction operator Here is the 
interaction strength between spins m and n. The coefficients 
{I3m,n] are treated as small quantities. In the absence of inter- 
action ({/3m,n) = 0), the propagator assumes a factorized form 

-^llA,,,,)=0 ~ -C|l)-C(2| ■ ■ ■ -L\N\, 

where X|„) = X(„| = Tr(p(„|e'°i"i) and 0(„| = aj'„. For 
{/0m,n) 0, we introduce additional factors (cluster correla- 
tions) to account for the interaction corrections. These cluster 
correlations can be introduced successively as follows. 

1. Two-spin correlations {-C)/,;)). 

If the bath consists of only two spins with indices {i, j], 
the propagator would be 

Xlui^Tr(p|,,|e'«''-^i) 

withpiyi =p(,) ®p|^] and 

Ov,,j\ = aiJj + ajP. +/3ij [JIJJ + jyj]) , 



i.e., X,{i,j} is obtained from Eq. (|4|l by dropping all spins 
except / and j. Without interaction (J3ij = 0), the prop- 
agator is 

The interaction correction makes the factorization to be 
- Thus the two-spin correlation is 

defined as 

X,,., - (7) 

Obviously, the Taylor expansion of the pair correlation 
with respect to the interaction strength is 

In = ciPij + C2/31 + ■ ■ ■ = 0(13), (8) 

where p denotes the typical magnitude of the interac- 
tion strength {yS,„,„). Thus In X.{i,j} is at most a first-order 
small quantity. 

2. Three-spin correlations {X.\ij^k\]- 

For a bath of three spins {;, j, k}, the propagator is 

XlU,«=Tr(p|y,„e"^"-«) 

with p|,j,i:) = p),) ®p|j| ®p^k] and 

m<« 

n=ij,k m,n=i,j,k 

i.e., X.{ij,k\ is obtained from Eq. Q by dropping all spins 
except i, j, and k. Similar to the two-spin case, 
can be factorized as 

where 



''' Ai\Aj]Ak]Ai.j]AjM-Cu,k) 

accounts for the non-factorizable correlation among the 
three spins. \nAi,j,k} would vanish if the interactions 
in <3(,j,A:) cannot connect the three spins {/, j, k} into a 
linked cluster. For example, if /?, ^ = jS, ^ = and /3j^k + 
0, the three spin propagator AuM would be factorized 
as 

AiJM = Ai}Aj,k\' 

which, together with In A'J] - 1" A'.k) - [according 
to Eq. dHJ], leads to InAijM - according to Eq. 
This connectivity property of lnX(i j,i:| leads to the Tay- 
lor expansion 

In Aij,k) = ciAvySa + c2/3j,if3j,k + CiPk^Pkj + Oiji"). (10) 
Thus In AiJ.k) is 3l most a second-order small quantity. 
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3. Cluster correlation {-Ccl- 

The above factorization procedure can be carried out for 
baths consisting of more and more spins. For a bath of 
an arbitrary group of n spins (denoted as C), the propa- 
gator becomes 



which is obtained from Eq. (|4| by dropping all spins 
except those belonging to the group C. By introducing 
the cluster correlation 



-Cc 



n -Cc 

C'cC 



the propagator is factorized as 
-Cc^Yllc'- 

C'CC 

By mathematical induction, it can be readily proved that 
In X.C vanishes if the interactions contained in Oq can- 
not connect all the spins in group C into a linked clus- 
ter. Such connectivity property ensures that in each 
term of the Taylor expansion of In about the inter- 
action strength, the coefficients {PijYs must appear at 
least (|C| - 1) times (|C| being the number of spins in the 
group). Thus InXc = O[0^^-^). 

In particular, the full propagator X of the whole bath is factor- 
ized into the product of all possible cluster correlations as 



CC(1,2,--,A'| 



An exact evaluation of the ensemble CCE in Eq. (fTTT l is 
not possible in general, which amounts to exactly solving the 
many-spin dynamics. In the qubit decoherence problem, it 
often suffices to truncate the CCE to an appropriate order M 
(denoted as M-CCE for short) by dropping all cluster correla- 
tions with sizes larger than M, 



(M) 



(12) 



\C\<M 



For example, the first-order truncation of the ensemble CCE 
(the 1-CCE) is 



\\£n 



(13) 



which is equivalent to Eq. (|4|l with all interaction terms in O 
dropped. In order to incorporate the interaction effects, the 
lowest nontrivial order of truncation is the second order (2- 
CCE), 



£ 



(2) 



i ' V {iji 



(14) 



which coincides with the pair-correlation approximationJ^i^S 
Since only connected clusters for which ln£c + con- 
tribute to the propagator £, the convergence (and hence the 
justification for the truncation) of the ensemble CCE can be 
estimated as follows. Fkst, if each spin interacts, on average, 
with q spins, then the number of connected size-M clusters is 
~ Nq'^^^ , with the total number of bath spins. Second, for 
a size-M cluster, InXc - 0(y6*^"'). The contribution to \n£ 

from all the size-M clusters is Yi\c\=m^^ £c ~ N{q/3)^^^. For 
q/3 < I, the ensemble CCE converges. 



B. General theory 

The above example can be readily generahzed to 

£^Tv{pe'^'e'^'-'---), (15) 

with a general ensemble p and an arbitrary series of time- 
ordered bath operators (9'-'* (j = 1, 2, • • ■ ). The bath interac- 
tions need not be purely off-diagonal or contain only pairwise 
interactions, and the initial density matrix need not be factor- 
izable. For a thermal ensemble as 



p = exp {-PrHj , 



(16) 



the density matrix itself can be viewed as a propagator with 
imaginary time r - -i/Sr and the whole propagator can be 
written as 



£^Tr{poe'("'^"h'^'\'^'''---), 



(17) 



with po denoting the trivial thermal state at infinite tempera- 
ture (jSr = 0). 

In essentially the same way as illustrated in the example 
above, a hierarchy of cluster correlations {Xc) can be intro- 
duced. First, the single-spin correlation is defined as 



£{i} = £{i). 



where 



X,„.Tr(p,„.<;.'^l?...) 



is obtained from Eq. ( fTSl l by dropping all spins except spin i 
and 



P[i} = Trk^i I 



(18) 



is the reduced density matrix of the /th spin. Then the cluster 
correlation for an arbitrary group C of bath spins is defined as 



£c^ 



£c 



n £c' 

C'cC 



where 



£c^Tr{pce'^ce'o'^' 
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is obtained from Eq. ( fTSl ) by dropping all bath spins outside 
group C and 



PC = Tiktc I 



(19) 



is the reduced density matrix of the cluster In particular, the 
whole bath propagator is factorized into all possible cluster 
correlations as 



x= n ^c, 

CC{\,2,--,N\ 



(20) 



which assumes exactly the same form as Eq. ( fTTT l. 

Now we discuss the convergence property of the ensemble 
CCE for general bath operators. To focus on the effect of the 
interaction strength, we consider the case that the initial bath 
density matrix is factorizable as in Eq. (|5]l. For a factoriz- 
able ensemble, it can be readily proved that lnX.c vanishes 
when the interactions contained in operators d^^\ d^^\ ■ ■ ■ 
cannot connect all the spins in the group C into a linked clus- 
ter. Considering that the operators (9'^^ contains the W-spin 
interaction such as B^^. . P- P for IV = 3 (in the example 

'l'2'3 '1 '2 '3 

of Sec. Ill Al W = 2), in the Taylor expansion of In X.c with 
respect to the dimensionless coupling coefficients yS's, the in- 



outside cluster C have been dropped. In contrast, the single- 
sample CCE in Ref. i24i keeps all the diagonal terms by re- 
placing the spins outside cluster C with their mean-field val- 
ues in the initial state For example, the diagonal in- 
teraction Y^nicJunJ^Jl between a spin J,ec inside cluster C 
and spins {intc) outside cluster C would be replaced with 
Jf Yjn<tc'yi,n(J'\J^,\J'}, which contributes a static local mean- 
field for the ith spin. With such a static mean-field treatment 
of the diagonal terms, the expansion is carried out with re- 
spect to the most essential dynamics of the spin bath, namely, 
the collective flip-flop of a cluster of spins that is responsible 
for the dynamical local field fluctuation for the qubit. This 
procedure, however, is not applicable to the ensemble CCE, 
for each sample state \ J'} from the ensemble p will generate a 
different static mean field. 

With the static local field fluctuation singled out, the clus- 
ter correlation in the single-sample CCE accounts for the 
collective dynamical local field fluctuation generated by off- 
diagonal interactions. As a result, its magnitude in single- 
sample CCE is determined by the magnitude of off-diagonal 
interactions, while in ensemble CCE it is determined by the 
magnitude of all kinds of bath interactions. As an example, 
consider pairwise interaction like 



teraction coefficients {B- 

connect all the spins in group C into a linked cluster. At least 
(M - 1) / {W - 1) interaction coefficients are needed to form 
a size-M linked cluster As a result, 



contained in each term must ^ = X Yj K'^" 



{M-\)l{W-\)\ 



(21) 



where yS is the typical magnitude of the coupling coefficients. 
The number of size-M clusters is ~ Nq^^\ So the total con- 
tribution to In X from all size-M clusters is 



\C\=M 

Therefore, the ensemble CCE converges for q^^^P < 1. In- 
terestingly, if the bath has certain initial correlations or entan- 
glement (such as a strongly correlated system at low tempera- 
ture), the convergence property would be determined by both 
the inverse temperature fSr and the typical coupling constant 



Let j6„d (J3d) be the typical magnitude of the off-diagonal (di- 
agonal) coupling coefficients j6"f„ (J3f„ „) and /3 be the greater 
one of /3d and j6nd- Then for a size-M cluster. In = Oifi'^^^) 
in the ensemble CCE, while InX^ = ^^nd^^^ single- 
sample CCE. The single-sample CCE would converge faster 
than the ensemble CCE. Moreover, the number of clusters in 
ensemble CCE is greater than that of the single-sample CCE. 
For example, if a spin J, in a cluster interacts with others 
through diagonal interaction only, then it generates no dynam- 
ical fluctuations and hence the single-sample cluster correla- 
tion vanishes, while the ensemble cluster correlation does not. 
For a bath with a relatively large number of spins and a factor- 
izable initial state, the ensemble CCE result would be close to 
the single-sample CCE using a random sampling of the initial 
bath state. The single-sample CCE is recommended in such 
cases. For a relatively small bath or a non-factorizable ensem- 
ble, the ensemble CCE is desirable. 



C. Comparison with single-sample CCE 

The evaluation of the single-sample average 

£^^(j\e^^'"e^^'"---\j) (22) 

on a factorizable product state \J'} - i8>„ \j„} can be done as 
a special case of the ensemble CCE by taking p = \J') {J'\ - 
®np\n] with = \i„){i„\. This method is slightly different 
from the previously developed single-sample CCE,^"* in defin- 
ing the cluster correlation for a cluster C. Here not only off- 
diagonal but also diagonal interaction terms involving spins 



III. NUMERICAL CHECK 

For a qubit-bath system described by a general pure dephas- 
ing Hamiltonian as in Eq. ([T]l, the decoherence of the qubit 
under the pulse control of the nth-order concatenated dynam- 
ical decoupling '^ii2i22i22i2£ is characterized by^° 

=Tr[pf/<->"'f/(+>], (23) 

where p is the density matrix for the initial bath state and tj^^'' 
are recursively defined as 
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Total evolution time T 



FIG. 1: (Color online) Qubit coherence in free-induction decay for 
a "sinusoidal" chain with A' = 500 spins: the exact solution (empty 
squares) vs. the results from ensemble CCE truncated to the first 
order (1-CCE, solid line) and the second order (2-CCE, dotted line). 

with u'^^ = e For example, free-induction decay, Hahn 

echo, and Carr-Purcell echo correspond to n = 0, 1, and 2, 
respectively. 

In this section, we consider an exactly solvable spin bath 
model (the one-dimensional spin- 1/2 XY model) and compare 
the qubit coherence X.„ from the ensemble CCE with the ex- 
act solutions by the Jordan- Wigner transformation.i^ii^ The 
A^-spin bath Hamiltonian conditioned on the qubit state |+) is 

H^-> ^±j]'^J;/flB„± y {JIJ- + , (24) 

«=i ^ «=i ^ ^ ' 

where z„ denotes the qubit-bath spin interaction strength (sim- 
ulating the hyperfine interaction strength for electron-nuclear 
spin systems), B„ is the intrinsic bath interaction strength, and 
bn is the interaction dependent on the qubit state. The bath is 
assumed to be in a high-temperature thermal ensemble with 
p = (1/2)'^. The qubit-bath interaction coefficients {z„} are 
taken from a sinusoidal distribution z„ = Zmax sin(nn/N) (re- 
ferred to as "sinusoidal" chain) or randomly selected from 
[0,Zmax] (referred to as "random" chain). Hereafter Zmax is 
taken as the unit of energy. The spin-flip interaction strengths 
{B„} and {b„} are randomly chosen from [10"-', 2 x 10"-'], cor- 
responding to typical bath spin flip-flop time r^f ~ 10^. The 
convergence of the ensemble CCE then requires B„T, b„T < 1 
or equivalently T < Tjf . 

First we consider the simplest case, namely, the qubit co- 
herence = Trj^pe'^* 'rg-/H'+'rj jn free-induction decay. The 
first-order truncation of the ensemble CCE gives 

X<'> = ncos(f), (25) 

i 

which is indeed the dephasing due to inhomogeneous broad- 
ening. In the short-time limit (z,T <K 1), Eq. ( |25] | becomes 




Pulse delay time T 




Pulse delay time T 




Pulse delay time T 




Pulse delay time T 



FIG. 2: (Color online) Qubit coherence in Hahn echo (n=l) and 
Carr-Purcell echo (n = 2) for different baths. The exact ensemble 
coherence Jl„ (empty squares) vs. the magnitude l-Cl'^'l of the exact 
single-sample coherence for three randomly chosen bath states \ J'), 
denoted by the solid, dashed, and dotted lines, respectively. 

^0) ^ g-r^rV2^ ^here 

is the variance of the "Overhauser" field h~ = Y^n^nJn of 
the spin bath. As evidenced by the good agreement between 
the 1-CCE and the exact result in Fig. [1] the ensemble free- 
induction decay is dominated by the inhomogeneous broad- 
ening, which leads to rapid decoherence within a time scale 
much shorter than the inverse qubit-bath interaction strength 
~ 1. The corrections due to spin-spin interactions, which 
could show up on a time scale comparable with the inverse 
interaction strength {l/B„,l/bn ~ 10^) and dominates the 
single-sample decoherence, is negligible during the ensemble 
free-induction decay process. 

To highlight the role of spin-spin interaction, we consider 
the qubit coherence X.„ (n - 1,2, ■ ■ ■) in spin echo or higher 
order concatenated control where the inhomogeneous broad- 
ening is eliminated and, consequently, the first-order trunca- 
tion of the ensemble CCE gives no decay: = 1. When 
the random phase factor leading to the inhomogeneous broad- 
ening is eliminated (as in the echo signals), the ensemble- 
averaged qubit coherence would be close to that averaged on 
a randomly sampled bath state, if the bath is relatively large. 
This is shown in Figs.|2](a) and|2](b) for a sinusoidal chain: 
The single-sample coherences for three randomly chosen bath 
states \ J') agree very well with the ensemble coherence. 

For a small spin bath, or for a random chain where the qubit 
decoherence is caused by the dynamics of a few small clus- 
ters, the qubit decoherence would depend sensitively on the 
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FIG. 3: (Color online) Ensemble qubit coherence in nth-order con- 
catenated dynamical decoupling for a sinusoidal chain with A' = 500 
spins: the exact ensemble solutions (empty squares) vs. the results 
from the ensemble CCE truncated to dilferent orders. 



choice of the initial bath state and hence the ensemble average 
would deviate significantly from the qubit coherence averaged 
on any specific sample of the initial state. This is clearly seen 
in Figs.|2](c) and|2](d) for a random chain. Even for a rel- 
atively large sinusoidal chain, the difference between single- 
sample decoherence and the ensemble average is noticeable 
when higher-order dynamical decoupling (e.g., Carr-Purcell 
echo) is applied [see Figs. |2je) and|2|f)]. This is because 
that under the higher-order control, the clusters responsible 
for the qubit decoherence grow larger and larger as the effects 
of smaller clusters are suppressed, and the specificity of the 
initial state of larger clusters is more important than that of 
smaller ones. 

Figure[3]compares the results from the ensemble CCE with 
the exact solutions for ensemble coherence under the control 
of concatenated dynamical decoupling of different orders. For 
the Hahn echo shown in Fig. [3ja), the second order trunca- 
tion of the ensemble CCE already agrees with the exact solu- 
tion very well, indicating that the decoherence is dominated 
by spin-pair dynamics. For Carr-Purcell echo in Fig. Ob), 
however, the second-order truncation becomes insufficient, as 



the leading order contributing from spin-pair dynamics has 
been eliminated and the coiTelated dynamics of larger clus- 
ters becomes important'^. Hence a fourth-order truncation of 
the ensemble CCE is required to reproduce the exact solu- 
tion. For concatenated dynamical decoupling of successively 
higher orders n, the leading contributions from successively 
larger clusters have been eliminated. As a result, successively 
higher-order truncations of the ensemble CCE are required to 
reproduce the exact results: 6-CCE for « = 3 in Fig. [Sjc) and 
8-CCE for « = 4 in Fig. Hd). 



IV. CONCLUSION 

The decoherence of a qubit in a spin bath is essentially 
determined by the many-body bath evolution, starting ei- 
ther from a pure or an ensemble state. As an extension 
of the previously developed single-sample cluster-correlation 
expansion^"* that addresses the spin bath dynamics starting 
from a noninteracting pure state, we have developed an en- 
semble CCE theory to solve the spin bath dynamics starting 
from an arbitrary ensemble state. In this approach, the en- 
semble propagator is factorized into the product of all possi- 
ble cluster correlations, each term accounting for the authen- 
tic (non-factorizable) collective excitation of a group of bath 
spins. For a finite-time evolution as in the qubit decoherence 
problem, convergent results can be obtained by truncating the 
ensemble CCE by keeping cluster correlations up to a certain 
size, as has been checked using an exactly solvable spin chain 
model (the one-dimensional spin- 1/2 XY model). For the con- 
vergence property (for a factorizable initial state), the ensem- 
ble CCE is determined by the typical strength of the bath spin 
interaction, being either diagonal or off-diagonal, while the 
single-sample CCE is determined only by the off-diagonal in- 
teraction coefficients. The ensemble CCE can be applied to 
baths with non-factorizable initial state as well. 
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